One-dimensional Rydberg Gas in a Magnetoelectric Trap 
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We study the quantum properties of Rydberg atoms in a magnetic Ioffe-Pritchard trap which is 
superimposed by a homogeneous electric field. Trapped Rydberg atoms can be created in long-lived 
electronic states exhibiting a permanent electric dipole moment of several hundred Debye. The 
resulting dipole-dipole interaction in conjunction with the radial confinement is demonstrated to 
give rise to an effectively one-dimensional ultracold Rydberg gas with a macroscopic interparticle 
distance. We derive analytical expressions for the electric dipole moment and the required linear 
density of Rydberg atoms. 

PACS numbers: 32.60.+i,32.10.Dk,34.20.Cf,32.80.Pj 



Because of their widely tunable properties, ultracold 
atomic gases provide the ideal playground to model and 
study complex many body systems. The interatomic in- 
teraction can be tailored using Feshbach-resonances and 
magnetic, optical, and electric fields can be applied in 
order to generate virtually any external potential. Fa- 
mous examples for the versatility are the demonstration 
of the Mott-Insulator to superfluid phase-transition of 
ultracold atoms in an optical lattice the BEC-BCS 
crossover in a gas of 6 Li or the Kosterlitz-Thouless 
phase transition studied within a two-dimensional Bose- 
Einstein condensate Traditionally, there is a great 
interest in studying systems with reduced spatial dimen- 
sions, such as the latter example. One paradigm is con- 
stituted by the work of Lieb and Liniger who were the 
first to solve the system of pointlike interacting bosons 
in one dimension using Bethe's ansatz In the limit 
of an infinitely strong interparticle interaction strength, 
a so-called Tonks- Girardeau gas emerges @. 

Besides gases of ground state atoms, particularly Ryd- 
berg gases represent excellent systems to study the influ- 
ence of a strong interparticle interaction on the dynamics 
of many-particle systems. Due to the large displacement 
of the ionic core and the valence electron, Rydberg atoms 
can develop a large electric dipole moment leading to a 
strong and long-ranged dipole-dipole interaction among 
them [6J. However, unlike ground state atoms Rydberg 
atoms suffer from radiative decay and hence the mutual 
interaction time is limited by the lifetime of the electron- 
ically excited state. But still, in so-called frozen Ryd- 
berg gases, where the timescale of the atomic motion is 
much longer than the radiative lifetime, exciting effects 
like the dipole-blockade Q or resonant population trans- 
fer [§| have been observed. 

Several works have focused on the important issue of 
trapping Rydberg atoms based on electric 0], optical 
f 101 ] , or strong magnetic fields ll( . Due to the high level 
density and the strong spectral fluctuations with spa- 



tially varying fields, trapping or manipulation in general 
is a delicate task. This holds particularly for the case 
where both the center of mass and internal motion are 
of quantum nature. Moreover, the inhomogeneous exter- 
nal fields lead to an inherent coupling of these motions. 
Addressing this regime, it has been recently theoreti- 
cally shown that Rydberg atoms can be tightly confined 
and prepared in long-lived electronic states in a magnetic 
Ioffe-Pritchard (IP) trap 1121] which can be miniaturized 
using so-called atomchips [13]. Here we use the IP con- 
figuration as a key ingredient in order to 'prepare' and 
study a one-dimensional (ID) Rydberg gas. Specifically, 
we propose a modified IP trap, a magneto-electric trap, 
which offers confining potential energy surfaces for the 
atomic center of mass (cm.) motion in which the atoms 
possess an oriented permanent electric dipole moment. 
We derive analytical expressions for the dipole-dipole in- 
teraction among trapped Rydberg atoms and estimate 
below which Rydberg atom density a ID Rydberg gas is 
expected to form. Moreover, we estimate the lifetime of 
such a gas. 



Proceeding along the lines of Refs. [12|, |14|, we em- 
ploy a two-body approach in order to model an alkali 
metal atom in a Rydberg state. We assume the single 
valence electron and the ionic core (mass M c ) to inter- 
act via a pure Coulomb potential. While the inclusion 
of the fine-structure and quantum defects can be read- 
ily done, it turns out not to be necessary for high an- 
gular momentum electronic states in the regime we are 
focussing on The IP field configuration is given by 
B(x) = Be 3 + Bi in (x) with Bi in (x) = G[x\ex - x 2 e 2 ] 
and the vector potential reads A(x) = A c (x) + Aii n (x) 
with A c (x) = §■ [x\e 2 — x 2 e\] and Ai in (x) = Gxix 2 e 3 , 
where B and G are the Ioffe field strength and the gra- 
dient, respectively. In addition, we apply a homogeneous 
electric field pointing in the xi-direction of the labora- 
tory frame F = Fe\. After introducing relative and 
cm. coordinates (r and R) and employing the unitary 



2 



transformation U = exp [i^ea x r ■ R] , the Hamiltonian 
describing the Rydberg atom becomes (atomic units are 
used unless stated otherwise) 



■ffiPE = H A + A c (r) p 



2M« 



/x JV -B(R) + F-r 



-/x e • B(R + r) + A lin (R + r) • p. 



(1) 



Here, Ha = *| — -is the Hamiltonian of a hydrogen atom 
possessing the energies E n = —^n~ 2 . The second term 
denotes the energy of the electron in the homogeneous 
Ioffe field due to its orbital motion. The following two 
terms of -f/rpE describe the motion of a point-like parti- 
cle possessing the magnetic moment fi N in the presence 
of the field B. The magnetic moments are connected 
to the electronic spin S and the nuclear spin S accord- 
ing to fi e — — S and fi N = — ^g-E, with g^ being the 
nuclear g-i actor. We neglect the term involving fj, N in 
the following due to the large nuclear mass. The electric 
field interaction, which in case of a neutral two-body sys- 
tem couples only to the relative coordinates, gives rise to 
the fifth term. The last two terms of Hipe are spin-field 
and motionally induced terms coupling the electronic and 
cm. dynamics. We focus on a parameter regime which 
allows us to neglect the diamagnetic interactions 

In order to find the stationary states of the Hamilto- 
nian ([T]), we assume that neither the magnetic nor the 
electric field causes couplings between electronic states 
with different principal quantum number n. In this case 
we can consider each n-manifold separately and may rep- 
resent the Hamiltonian ([!} in the space of the 2n 2 states 
which span the n-manifold under investigation. The pa- 
rameter range in which this approximation is valid has 
been thoroughly discussed in Refs. II, 15 1. Because of 
the translational symmetry of the IP and the electric 
field, the axial cm. motion along Z can be separated from 
the transversal motion in the X-Y plane. If we omit the 
energy offset E n and introduce scaled cm. coordinates 
(R — * 7~aR with 7 = GM C ) while scaling the energy 
unit with e sca i e = 73 /M c , we arrive at the Hamiltonian 



H 



Pi 



H e 



(2) 



This Hamiltonian governs the transversal cm. as well 
as the electronic dynamics and involves the effective 
magnetic field G = Xe x - Ye 2 + 7"lM c Be 3 . The 
symbols fi, H m , and H e are the 2n 2 -dimensional ma- 
trix representations of the operators \ [L + 2S], H m = 
An n (r) • p + Bij n (r) • S, and the electric field interaction 
H e = Fx, respectively (we introduced L = r x p). A 
thorough interpretation of this Hamiltonian in the ab- 
sence of H e is provided in Ref. 12]. In order to solve 



the corresponding Schrodinger equation we employ an 
adiabatic approach. To this end an unitary transforma- 
tion U(X, Y) which diagonalizes the last two (matrix) 
terms of the Hamiltonian is applied, W(X 7 Y)(fi ■ G + 



6.5 



6.4 



6.3 



6.2 



6.1 



o 6 



X 

s 

LU 



(a) 



89.9 MHz 




x (b) A 

0.6 MHz 



v(e) 



(f) 



6.5195 
6.5185 

6.3825 
6.382 



6.301 
6.3005 

6.1695 
6.169 



-30000 30000 x 
X (a.u.) 

FIG. 1: (Color online) Potential energy surfaces of the cm. 
motion of a 8r Rb atom [n = 30) in an IP trap with B = 
10 G, G = 10 Tm" 1 . Dashed lines: F = 0, solid lines: F = 
5.14Vm-\ An overview of the seven energetically highest 
potential curves is shown in panel (a) . Magnified views of the 
uppermost (b,c) and next lower ones (d,e) are also provided. 
The range of the X-coordinate, corresponding to 2.1/^m, is the 
same for each subfigure (a)-(e). The total field configuration 
is sketched in panel (f ) where the circles depict the locations of 
the minima of the uppermost (big circle) and the two adjacent 
lower-lying (small circles) cm. surfaces. The magnetic field 
lines are indicated in gray while the electric field is sketched 
by black arrows. 



j-3M c [H m + H e ])U(X,Y) = E a {X,Y). Since U(X,Y) 
depends on the cm. coordinates, the transformed ki- 
netic energy term involves non-adiabatic (off-diagonal) 
coupling terms which can be neglected in our parameter 
regime 
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We are thereby led to a set of 2n 2 decoupled 
differential equations governing the adiabatic cm. mo- 
tion within the individual two-dimensional energy sur- 
faces E a (X,Y), i.e., the surfaces E a (X,Y) serve as po- 
tentials for the cm. motion of the atom. In Fig. [1] we 
present intersections along the X-direction of such poten- 
tial surfaces for B = 10 G, G = lOTrn -1 and n = 30 in 
the case of 87 Rb. For zero electric field strength (dashed 
lines) , the potential curves are organized in groups which 
are energetically well-separated by a gap of y~sM c B = 
89.9 MHz. The uppermost surface is non-degenerate and 
provides an approximately harmonic confinement with a 
trap frequency of u> = Gy/n/2BM c — 13.9 kHz corre- 
sponding to 0.1 /iK. The two adjacent lower surfaces are 
degenerate and also approximately harmonic. As soon 
as an electric field is applied, all surfaces are shifted 
considerably in energy. This is visible from the solid 
curves in Fig. Q] for which an electric field of strength 
F = 5.14 Vm _1 is applied. The shapes of the potentials 
are barely affected by the electric field such that Rydberg 



states which were trapped in a pure IP configuration re- 
main confined also in the magneto-electric trap. More- 
over, adding the electric field leads to non-trivial effects: 
The second and third surface, which were almost degen- 
erate in the absence of the electric field, are now shifted in 
opposite ways along the xi-direction. All surfaces shown 
provide a harmonic confinement with a trap frequency 
u> also in the ^-direction. We remark that the chosen 
parameter set does not generate an extreme constella- 
tion hence an even stronger confinement can be achieved 
without invalidating the applied approximations (l^ j. 

Let us now investigate the electronic properties of a 
Rydberg atom being trapped in the uppermost potential 
surface. For F = and sufficiently large values of B, this 
surface is formed almost exclusively by the highest pos- 
sible electronic angular momentum state, i.e., I — n — 1 
(l2T |. If F is increased, electronic states with smaller I will 
be inevitably admixed to the electronic state belonging 
to this energy surface. An interesting property to in- 
vestigate is hence the electric dipole moment of trapped 
Rydberg atoms: While for F — the electronic states 
are almost pure parity eigenstates and therefore exhibit 
almost no electric dipole, the admixture of lower I states 
to the uppermost surface in the presence of the field is 
expected to give rise to a non- vanishing expectation value 
of the dipole operator. Indeed, this becomes evident in 
Fig. [5] where the uppermost potential surface and the 
three components of the expectation value of the elec- 
tric dipole operator D(R) = (r) (R) are shown (same 
parameters as in Fig. It can be clearly seen that a 
permanent dipole moment is established whose dominant 
contribution points along the electric field vector. 

In order to study the dependence of D(R) on the field 
strengths F and B as well as on the degree of elec- 
tronic excitation, we use perturbation theory. In the 
limit of a large Ioffe field strength B the unitary trans- 
formation which diagonalizes the Hamiltonian J5J can be 
written explicitly as U T = e -i a (L*+S x ) e -ip{L v +s v ) 
sina = -r|G|-\ cosa = y/jGf^Y 2 |G|"\ sin/3 = 
X(\G\ 2 -Y 2 )-i and cos/3 = >y-*M c B (|G| 2 - Y 2 )~^ . 
This transformation rotates the z-axis into the local mag- 
netic field direction where a and f3 denote the rotation 
angles. Using this result we find up to first order in F/B 
the electric dipole moment (in atomic units) 



D(R) 
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n 2 (n — 1) 



cos/3 
sin P sin a 
sin (3 cos a 



(3) 



We note that D(R) scales proportional to the third 
power of the principal quantum number and can there- 
fore gain a significant magnitude even if the ratio F/B 
is small. Good agreement of Eq. ([3|) with the calculated 
data presented in Fig. [2] is found; e.g., in the vicinity of 
the minimum of the potential surface (X = Y = 0) we 
find an exact value of D x — 270 whereas the expression 
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FIG. 2: (Color online) (a) Uppermost electronic potential sur- 
face for the cm. motion of 87 Rb in the n = 30 multiplet 
and the parameters used in Fig. [T] (b-d) Components of the 
electronic dipole moment D(R) in atomic units. One recog- 
nizes the clear alignment of the electric dipole moment along 
the electric field vector. The numerically calculated values of 
D(R) are to good accuracy reproduced by Eq. ©. 



© yields 276. For the remaining components Eq. ([3]) 
yields zero at the origin. For smaller ratios of F/B, even 
better agreement can be achieved. 

Because of the dependence on the angles a and /3, the 
dipole moment depends weakly on the quantum state of 
the cm. motion. However, since the field configuration 
is translationally symmetric the electric dipole moment 
is independent of the Z-position of the Rydberg atoms in 
the trap. If we now consider two trans versally confined 
atoms in the same trap at the longitudinal positions 
and Zg , we can write for their dipole-dipole interaction 

U D (R.4, R B ) = — -L— 3 (D(R A ) • D(R B ) 

-3[D(R A ).e][D(R B ).e])«H(M^) (4) 

where e denotes the interparticle unit vector. The ap- 
proximation in Eq. (U) holds due to the orientation of 
the dipoles and the assumption that \Za — Zb\ is large 
compared to the transversal oscillator length of the trap. 
These conditions moreover ensure a minimal coupling of 
the transversal and longitudinal motion. Using this ap- 
proximation one can estimate the interaction energy of 
one atom being part of an infinite atomic chain with an 
interparticle spacing a. One finds 



2 2!<2>f>- 



k=l 



81 F 2 
2a3~B 2 



n 4 (n-I) 2 C(3) (5) 



with the Ricmann zeta function £(x). Here we have ap- 
proximated D 2 (R) as D 2 (0) since the dipole moment 
barely varies in the vicinity of X = Y = 0. If the in- 
teraction energy E- ln t is smaller than the transversal trap 
frequency ui, we can assume that the interacting atoms 
remain in the transversal ground state: This is consid- 
ered the ID regime. The linear density below which a ID 
Rydberg gas is expected to form is then given by 
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C(3)n y / 2 (n-I) 2 



(6) 
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Above this density, excited transversal cm. states might 
be populated resulting in a quasi ID Rydberg gas which 
is certainly of interest on its own. For our parameter set, 
we obtain a minimal interparticle spacing a = 43 /im; 
hence a chain of 1 mm in length contains 23 particles. 
This density can be further increased by either increasing 
the magnetic field gradient and/or decreasing the electric 
field strength: At B = 10 G, G = IOOTitT 1 , and F = 
0.514 Vm _1 a chain of the same length would contain 230 
Rydberg atoms. 

Finally, the issue of radiative decay has to be ad- 
dressed. Since the electric field admixes merely a few 
I = m states (I < n — 1) to the electronic wave func- 
tion of the uppermost surface, its circular character re- 
mains dominant resulting in one prevalent decay chan- 
nel. For an atom being confined to the energy surface 
which is shown in Fig. [51 we have calculated a lifetime 
of r w 2.1ms which is in good agreement with the cor- 
responding field-free result r(n,n — 1) « ^ (a) 5 Hi- 
Corrections to this bare decay rate are found to be of the 
order of (F/B) 2 n 3 . Due to the scaling proportional to 
n 5 , the lifetime can be significantly enhanced by exciting 
to a higher principal quantum number n. In addition, 
it can be further prolonged by establishing an adapted 
experimental setup which inhibits the electromagnetic 
field mode at the dominant transition frequency 17 1. 



At the same time, a cryogenic environment will dimin- 
ish the undesirable effect of stimulated (de-)excitation 
by blackbody radiation. The timescale of the dynam- 
ics of the Rydberg chain on the other hand depends on 
the field strengths via the dipole moment and the in- 
terparticle spacing: A harmonic approximation of the 
dipole-dipole interaction yields the one-particle oscilla- 
tor frequency uj^d = \/24 D 2 (0)/M c a 5 . As an example, 
the field configuration B = 10 G, G = lOOTm -1 , and 
F = 0.514 Vm -1 yields a timescale of less than 1 ms. 

Let us now briefly comment on the realization of such a 
Rydberg gas which is certainly a challenging experimen- 
tal task. One could start from an extremely dilute ultra- 
cold atomic gas prepared in an elongated IP trap. For 
transferring ground state atoms to high angular momen- 
tum Rydberg states, techniques such as crossed electric 
and magnetic fields or rotating microwave fields can be 
employed, see Ref. [18J and references therein. For low 
angular momentum states, trapping and the formation 
of a permanent dipole represent still an open question 
since quantum defects, spin-orbit coupling and reduced 
radiative lifetimes have to be taken into account. During 
the preparation, the excitation lasers have to be focussed 
such that Rydberg atoms emerge only at positions sepa- 
rated by the interparticle spacing a which is required to 
meet the criterion ([6]) . Since a is in the order of several 
/xm, which can be resolved optically, this should be feasi- 
ble. The large value of a moreover ensures that the mu- 
tual ionization due to the overlap of the electronic clouds 
of two atoms does not occur. For our circular states 



with n = 30, the atomic extension can be estimated by 
(r) « n 2 = 48 nm and is thus orders of magnitude smaller 
than the corresponding value of a for our field configu- 
ration. In order to probe the dynamics of the resultant 
Rydberg chain, one can field-ionize the atoms: From the 
spatially resolved electron signal a direct mapping to the 
positions of the Rydberg atoms should be possible. 

In conclusion, we demonstrated that in a magneto- 
electric trap Rydberg states can be confined in electronic 
states exhibiting a permanent electric dipole moment of 
hundreds of Dcbyes. Analytical expressions for the den- 
sity which is required to enter the ID regime were calcu- 
lated. Moreover, we pointed out that the lifetime of the 
Rydberg states is sufficiently long to probe the dynamics 
of the interacting gas. This regime is complementary to 
the well-studied frozen Rydberg gases where mechanical 
atom-atom interaction effects can hardly be probed. The 
potential of the proposed magneto-electric trap is by far 
not entirely exhausted; e.g., one could think of using the 
double- well structure which is visible in Fig. QJd) m order 
to realize two coupled dipolar Rydberg chains. 
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